Pygmy dipole resonance in ^^'Pb 
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Scattering of protons of several hundred MeV is a promising new spectroscopic tool for the study 
of electric dipole strength in nuclei. A case study of ^"^Pb shows that at very forward angles J'^ = 1~ 
states are strongly populated via Coulomb excitation. A separation from nuclear excitation of other 
modes is achieved by a multipole decomposition analysis of the experimental cross sections based 
on theoretical angular distributions calculated within the quasiparticle-phonon model. The B(i?l) 
transition strength distribution is extracted for excitation energies up to 9 MeV, i.e., in the region 
of the so-called pygmy dipole resonance (FDR). The Coulomb-nuclear interference shows sensitivity 
to the underlying structure of the El transitions, which allows for the first time an experimental 
extraction of the electromagnetic transition strength and the energy centroid of the FDR. 

PACS numbers: 25.40.Ep, 21.10.Re, 21.60.Jz, 27.80.-fw 



The electric dipole (-El) strength in nuclei is dominated 
by the isovector giant dipole resonance (GDR), originat- 
ing from the collective motion of neutrons against pro- 
tons in the nucleus and located well above the particle 
emission threshold. The GDR provides basic insight into 
the isovector properties of the nuclear force and thus was 
intensively investigated both experimentally and theo- 
retically [ll, 1^. However, at present the interest is more 
focused on low-lying dipole strength, well below the GDR 
energies, referred to as pygmy dipole resonance. R ap- 
pears in nuclei with neutron excess and might be pictured 
macroscopically to result from oscillations of these excess 
neutrons against an inert core with N ^ Z (see, e.g., 
Ref. [3| and refs. therein). The PDR is predicted in all 
microscopic calculations based on the random-phase ap- 
proximation (RPA) , but the predicted central energy and 
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strength often differ considerably, in particular between 
those based on non-relativistic and relativistic mean-field 
approaches (see e.g. Ref. [4] for an example in the stable 
tin isotopes). 

By its nature, the PDR may shed light onto the forma- 
tion of neutron skins in nuclei [5|-l7| and, because of the 
strong correlation in the RPA models, in turn on the sym- 
metry energy [a, |8|, |9| . Constraints on the magnitude and 
density dependence of the symmetry energy are impor- 
tant ingredients for the modeling of neutron stars. There 
is a clear correlation between the total electric dipole po- 
larizability and the neutron skin [9, jL(|. However, it has 
been argued that the PDR alone carries independent in- 
formation [q, '95 . 

The properties of the PDR in stable nuclei have been 
studied extensively for different neutron and proton shell 
closures with the (7, 7') reaction (e.g. Ref. [11] and ref- 
erences therein). Crucial data in exotic neutron-rich nu- 
clei, where the PDR should be enhanced, are still scarce 
BllJ, ll3l ^"^d suffer from large systematic uncertainties. 



The heaviest stable doubly magic nucleus ^"^Pb has al- 
ways been a benchmark and the PDR has been studied 
in various recent (7,7') experiments [ij-[l4|- Theoreti- 
cally, closed-shell nuclei permit the inclusion of complex 
degrees of freedom beyond the mean-field level in the mi- 
croscopic calculations [ij, [l^ . 

While the results of Refs. |14l - [l7| agree quite well, 
the (7,7') reaction in general suffers from two problems: 
the experimental quantity measured is the ground-state 
gamma decay width times the ground-state branching ra- 
tio. The latter is usually not known and assumed to be 
100% in most analysis. However, statistical model cal- 
culations indicate potentially large correction factors (al- 
though not for the case of ^°*Pb) modifying the resulting 
PDR strength distributions considerably [19J. Further- 
more, the dominance of particle over 7 decay suppresses 
the experimental signal above threshold. Another, more 
general problem is the experimental separation of the 
PDR from the GDR. While theoretical transition den- 
sities provide a signature in the model calculations |20| . 
the experimentally determined reduced B(i51) strengths 
do not allow such a distinction. A possible experimental 
approach to the distinction of PDR and GDR is the use 
of isoscalar probes, and a pioneering (a, a'7) experiment 
has been performed [2l| demonstrating a large exhaus- 
tion of the isoscalar energy-weighted sum rule by low- 
energy El transitions in ^"^^Pb. 

Recently, a new experimental technique utilizing po- 
larized proton scattering at and close to 0° to measure 
the complete El strength in nuclei has been developed 
[22l |. It allows, in particular, a consistent extraction of 
the El transition strengths below and above the neu- 
tron threshold. At small momentum transfers and inci- 
dent proton energies of several hundred MeV the cross 
sections are dominated by isovector spinflip-Ml transi- 
tions (the analog of the Gamow- Teller mode) 23] and 
by Coulomb excitation of non-spinflip El transitions. A 
separation of these two contributions can be achieved ei- 
ther by a multipole decomposition analysis (MDA) of the 
angular distributions or by the analysis of polarization 
transfer observables. In Ref. [2J] excellent agreement of 
the two methods was demonstrated for the case of ^°^Pb. 
Here we present our results for the El strength distri- 
bution at energies below the GDR based on the MDA. 
We also demonstrate that because of the interference of 
Coulomb and nuclear interaction, the angular distribu- 
tions do show sensitivity to the underlying structure al- 
lowing for an experimental separation of PDR and GDR 
contributions. 

The ^'^^Pb(p,p') experiment was performed at the 
RING cyclotron facility of the Research Center for Nu- 
clear Physics (RCNP), Osaka University. A description 
of the experimental technique can be found in Ref. [22| 
and details of the present experiment and the polariza- 
tion transfer measurements in Ref. [2^. A proton beam 
of 295 MeV with intensities 2 — 10 nA and an average 
polarization of Pq ~ 0.7 bombarded an isotopically en- 
riched ^°^Pb foil with an areal density of 5.2 mg/cm^. 



Data were taken with the Grand Raiden spectrometer 
[26l | in an angular range 0° — 2.5° and for excitation en- 
ergies i^a; ~ 4 — 22 MeV. Additional data with unpolar- 
ized protons were taken at angles up to 10°. Employ- 
ing dispersion matching techniques, an energy resolution 
AE ~ 25 keV (full width at half maximum) could be 
achieved. A spectrum of the ^°^Pb(p,p') reaction with 
the spectrometer set at 0° is shown in Fig. [1] for the ex- 
citation energy range 4.5 — 9 MeV, where the PDR is 
expected to lie. The arrows indicate the excitation en- 
ergy of excited states identified in the (7,7') experiments 
14l-[l7[. Essentially, all prominent dipole transitions ob- 
served in the latter experiments are also excited in the 
present measurements. 
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FIG. 1: (Color online) Low-energy part of the spectrum of the 
^"*Pb(p,p') reaction at Ep = 295 MeV and Qiab = 0° - 2.5°. 
The arrows indicate transitions also observed in ^''*Pb(7,7') 



experiments [14H1 

Excitation of 1~ states is possible through nuclear and 
Coulomb interaction, and both contributions add coher- 
ently to the cross sections. To verify the assumption of 
a predominant Coulomb excitation at angles close to 0°, 
predictions for the angular distributions in {p,p') scat- 
tering were calculated based on a semiclassical model 
[27| . As examples, results for the prominent transitions 
to 1" states at E^ = 5.512 MeV and 6.720 MeV are 
shown in Fig. [21 Because of the finite angular resolution 
of the Grand Raiden spectrometer, the calculated cross 
section angular distributions were convoluted with Gaus- 
sian functions with widths corresponding to the vertical 
and horizontal angular opening of the detector system. 
The shape of the experimental angular distributions is 
well described and their absolute magnitudes can be re- 
produced when the calculations are normalized to the 
average B(E1) strengths deduced from the (7,7') experi- 
ments [14l4l7l | . The remaining deviations at angles larger 
than 2° are attributed to effects of Coulomb-nuclear in- 
terference and in case of the transition to the state at 
6.720 MeV to contributions from unresolved transitions 



with higher multipolarities. 
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FIG. 2: (Color online) Angular distributions of the promi- 
nently excited 1" states at Ex = 5.512 MeV and 6.719 MeV in 
the '^°^Ph{p,p') reaction at Ep = 295 MeV. The dashed lines 
are predictions of Coulomb excitation cross sections based on 
the semiclassical approach [27l |. 

In order to determine such contributions and to en- 
able a separation from the smn-Ml resonance known 
to set in at E^ > 7 MeV [ll, [111, a MDA was per- 
formed. The method, based on model predictions of 
the angular distribution shapes, is commonly used in the 
analysis of complex spectra from hadronic reactions, e.g. 
for an extraction of B(GT) strengths in charge-exchange 
reactions |29| or isoscalar giant resonance strength dis- 
tributions from inelastic a-particle scattering [30|, and 
also for inelastic electron scattering form factors of nu- 



clei [3l|. Theoretical proton scattering cross sections 



were calculated using the code DWBA07 |32| with RPA 
amplitudes and single-particle wave functions from the 
quasiparticle-phonon model (QPM) [JM . The t- matrix 
parametrization of Love and Franey [33| at 325 MeV was 
used as effective projectile-target interaction. For each 
discrete transition {E^ < 7 MeV) or excitation energy 
bin {Ex > 7 MeV, cf. Tab. |l|, the experimentally ob- 
tained angular distributions were fitted by means of the 
least-square method to a sum of the calculated angular 
distributions weighted with coefficients a^/^^ (with the 
condition a^/*'^^ > 0). 

Some approximations were necessary to make the 
MDA tractable: Experimental data, although available 
up to 10°, were restricted to scattering angles Qiab < 4° 
because of the increasing complexity of contributions 
from different multipoles at higher momentum transfers 
34| . Isovector spin-Ml excitations were represented by 
a single characteristic curve, justified by the identical an- 
gular dependence of the cross section for all transitions 
of this type in the angular range considered. Further- 
more, only El transitions with a strength larger than 
0.01 e^fm^ were taken into account. Figure EJa) com- 
pares the shape of isovector spin-Af 1 transitions with 



representative examples of E\ transitions to states of 
the FDR and GDR, respectively. Indeed, the latter can 
not only be distinguished from the Ml case but also be- 
tween each other. All other contributions to the cross 
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FIG. 3: (Color online) Comparison of theoretical angular dis- 
tribution shapes used in the MDA: (a) Ml and representative 
El transitions to the FDR and GDR, (b) E2 and EA, (c) EZ, 
M2 and A/3. All curves are normalized at 0iab ~ 1.5°. 



sections were substituted by angular distributions of ei- 
ther E2 or £'3 transitions, whose shapes were taken to 
be that of the most collective transition of each type. 
Other multipolarities of potential relevance like Af2, A'/3 
or £'4 exhibit very similar angular distribution shapes to 
either E2 or E3 as shown in Figs. Eljb) and (c). Since 
isoscalar monopole transitions are only weakly excited in 
proton scattering and the giant monopole resonance is 
located at higher excitation energies than the region of 
interest, possible contributions from EO transitions were 
neglected. 

The final coefficients were obtained by computing the 
MDA for all possible combinations of El, Ml and E2 
(or E3) transitions and taking the x^-weighted average 
of all individual a^/^^'*' values. Examples of fits for two 
adjacent energy bins in the region of overlapping levels 
[see Figs. Hl^a) and (b)] demonstrate the sensitivity of 
the MDA to distinguish El and A/1 contributions to the 
cross sections with a dominance of El in (a) and of A/1 
in (b), respectively. 

Table U summarizes the MDA results for excitation en- 
ergies Ex =4.8-9 MeV. The partial El and A/1 cross 
sections listed are integrated over a scattering angle range 
0° — 0.94°, where nuclear contributions to the El cross 
sections are on the level of 1% only. Thus, B(/?l) tran- 
sition strengths can be directly derived. 

As discussed in Ref. [ij|, the QPM calculations used 
for the present MDA analysis can distinguish between 
El transitions to the PDR and GDR by an analysis 
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TABLE I: El and spin-A/1 cross sections in the ^°*Pb(p,p') 
reaction for excitation energies Ex ~ 4.8 — 9 MeV integrated 
over scattering angles 0° —0.94° and derived B(_E1) strengths. 



FIG. 4: (Color online) Examples of the MDA fits for two 
adjacent energy bins in the energy region of overlapping levels. 



of the theoretical transition densities. The examples in 
Fig. EJa) demonstrate that the resulting {p,p') angular 
distributions also show significant differences between ex- 
citation of these two modes. Thus, the MDA can pro- 
vide information on the dominant structure of the El 
transitions. In Fig. [5l the best x^ values of least-square 
fits using either PDR (blue diamonds) or GDR (red di- 
amonds) transition amplitudes are plotted as a function 
of the excitation energy. For E^ < 8.23 MeV, x^ i'^- 
sults assuming a PDR structure of the El excitations 
are consistently superior to fits with GDR-type angu- 
lar distributions, while the reverse is observed for higher 
excitation energies. This finding allows to experimen- 
tally extract the properties of the PDR in ^"^Pb, viz. 
a centroid energy E ^mc — 6.82(2) MeV and a strength 
X;B(£'1) = 2.18(7) e^fm^. The systematic uncertainty 
of the MDA approach is estimated to be at most 10% 
by constructing mixed PDR/GDR angular distributions 
with given amplitude ratios and studying the variation 
of x^ in the fits. 
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FIG. 5: (Color online) Best x^ values in the MDA using ei- 
ther PDR- or GDR-type angular distributions for excitation 
energies E^ ~1 — 'i MeV and bins defined in Tab. HI 

Figure|5]compares the electric dipole strength distribu- 



tion up to i^a; = 9 MeV extracted from the present exper- 
iment (a) with previous results combining (7,7') [14l4l7l | 
and 207pb(n,7)_Ip data up to 8 MeV with a 208pb(e,e') 
measurement |36j at higher excitation energies (b). The 
agreement is excellent up to the neutron separation en- 
ergy (Sn — 7.33 MeV). Above threshold the present work 
shows additional, previously unobserved strength. This 
can most likely be attributed either to unknown neutron 
partial decay widths of the excited states. Above 8 MeV, 
El strength from the present work is consistent with the 
results of Ref. [3g| within the systematic uncertainties of 
the latter. 

Comparison with theoretical results from the 
quasiparticle-phonon model (QPM) [iJI, the rela- 
tivistic time-blocking approximation (RTBA) [37|, and 
a shell-model (SM) calculation based on the interaction 
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FIG. 6: (Color online) El strength distributions in ^°*Pb be- 
tween 4.8 MeV and 9 MeV from (a) the present experiment 
in comparison with (b) previous results and theoretical calcu- 
lations within the (c) QPM, (d) RTBA and (e) sheU model. 
Note the scale reduction by a factor of 10 for the RTBA re- 
sults. 



described in Ref. [38] is displayed in Figs. M,c)-{e), re- 
spectively. The QPM and RTBA approaches have been 
discussed in detail in Ref. [24] and the SM in Ref. fl7l |. 
The QPM, including a model space up to 3-phonon 
states, provides a good description up to 7.5 MeV. At 
higher excitation energies the strength is clearly too 
small. The RTBA calculation so far includes beyond 
the 1 particle- 1 hole {Iplh) states only configurations 
of the type Iplh ® phonon, and thus the fragmentation 
is insufficient to describe the data (note the general 
reduction by a factor of 10 in Fig. [BUd) to bring it on 
the scale of the other results). The total strength up 
to 9 MeV is more than a factor two too large. The SM 
shows a fair agreement for the total strength but the 
fragmentation is somewhat underpredicted. 



It would be interesting to extract from the models cen- 
troid and strength of the PDR in ^°*Pb to be compared 
with the experimental result quoted above. However, 
we refrain from giving the corresponding values summed 
over the interval where the PDR is found experimentally. 
In the models the properties of the PDR are very sensi- 
tive to the mean-field description and predicted strength 
might thus be partially shifted to higher excitation ener- 
gies. One should rather analyze the theoretical transition 
densities [20, [33| to select those with a dominant PDR 
character. 

To summarize, high-energy proton scattering at angles 
close to and including 0° is used as a new experimental 
method to determine the B(£'l) strength distribution be- 
low the GDR in ^'^^Pb. Combined with dispersion match- 
ing techniques, it provides a novel spectroscopic tool for 
the study of the electric dipole response in nuclei. The 
method also overcomes some limitations of other com- 
monly used experimental techniques like (7,7') or (7,^) 
reactions, which are sensitive to assumptions about decay 
branching ratios and typically limited to energy regions 
either below or above the neutron separation energy. 

The MDA exhibits sensitivity to the structure of El 
transitions in ^°*Pb through the difference of Coulomb- 
nuclear interference contributions to the cross sections of 
PDR and GDR excitations. This is an important finding 
because the strength distribution does not allow a dis- 
tinction between the different modes. Thereby, the elec- 
tromagnetic transition strength and the centroid energy 
of the PDR could be determined for the first time, pro- 
viding an experimental benchmark for the strongly vary- 
ing theoretical predictions (cf. Fig. |6]). The method is 
not restricted to closed-shell nuclei but can be applied in 
vibrational and well deformed nuclei where microscopic 
calculations provide a good description of nuclear exci- 
tations. Possible approaches how to extract the PDR 
content from model calculations of the El response have 
been discussed e.g. in Refs. [40, |41|. Alternative infor- 
mation on the PDR structure may be obtained from 
isoscalar probes 42, 43] as demonstrated by studies of 
the (a, a'7) reaction |2l U 44 | . 
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